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Abstract 


The main objective of this paper is to develop a seven-link dynamic model of the operator’s body while 
working with a motorized backpack sprayer. This model includes the coordinates of the sprayer relative to the 
body, the rotational inertia of the sprayer, the muscle moments acting on the joints, and a kinematic coupling that 
keeps the body balanced between the two legs. The constraint functions were determined and the non-linear 
differential equations of motion were derived using Lagrangian equations. The results show that undesirable 
fluctuations in the ankle force are noticeable at the beginning and end of a swing phase. Therefore, injuries to the 
ankle joint are more likely due to vibrations. The effects of engine speed and sprayer mass on the hip and ankle 
joint forces were then investigated. It is found that the engine speed and sprayer mass have significant effects on 
the hip and ankle forces and can be used as effective control parameters. The results of the analysis also show 
that increasing the engine speed increases the frequency of the hip joint force. However, no significant effects on 
the frequency of the ankle joint force are observed. The results of this study may provide researchers with insight 
into estimating the allowable working hours with the motorized backpack sprayers, prosthesis design, and load 
calculations of hip implants in the future. 


Keywords: Lagrange equation, Operator, Vibration, Weight of sprayer 


Introduction 


One of the most popular ways of crop 
protection against weeds and pests is applying 
motorized backpack sprayers. The use of the 
motorized backpack sprayer eliminates the 
need for hand pumping and is suitable for 
small-scale farms. However, the major 
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disadvantage of motorized backpack sprayers 
is the external forces acting on the operator’s 
body (Kouchakzadeh, & Beigzadeh, 2015). 

Forces acting on the human body are 
important factors in the initiation and 
progression of joint disease (Astephen, 
Deluzio, Caldwell, & Dunbar, 2008). Force 
analysis of hip and ankle joints can be useful 
in the development of strategies to avoid and 
manage conditions such as osteoarthritis and 
deterioration of femoroacetabular (Correa, 
Crossley, Kim, & Pandy, 2010). Force analysis 
of hip and ankle joints requires multi-segment 
models. 


2 Journal of Agricultural Machinery Vol. 14, No. 1, Spring 2024 


Several researches have been done on the 
dynamic modeling of the human body. Kuo 
(2001) applied a simple model of bipedal 
walking to evaluate simple hypotheses for the 
metabolic cost of muscle activity. Tlalolini, 
Chevallereau, & Aoustin (2010) optimized the 
walking motions of a bipedal model by 
minimizing the value of the torque. The 
optimization process was carried out with and 
without the rotation of the supporting foot 
about the toe. Huang, Wang, Chen, Xie, & 
Wang (2012) employed a seven-link dynamic 
bipedal walking model with flat feet to analyze 
the dynamics of walking. The results indicated 
that ankle stiffness plays different roles in 
different gaits. Martin & Schmiedeler (2014) 
proposed four and six link planar biped models 
with knees and rigid circular feet; The ankle 
was not considered in the four link model. The 
results showed that the model with ankles is 
more accurate in predicting energy 
consumption during normal human walking at 
different speeds. Sharbafi & Seyfarth (2015) 
utilized a five link model with a rigid upper 
body and two segmented legs to extract 
internal relations between the joints’ angles 
and stance leg orientation which hold the 
configuration harmonized during the gaits. 
Jena, Kumar, Singh, & Mani (2016) developed 
a biomechanical model to predict metabolic 
energy consumption for carrying load 
manually by varying modes (head, shoulder, 
and back), loads, and ground inclinations. The 
results indicated that frontal torque (in 
shoulder mode) requires higher physiological 
energy than sagittal torque. Ma, Xu, Fang, Lv, 
& Zhang (2022) present the dynamic model of 
the human-prosthesis heterogeneous coupled 
system. Ma, Xu, & Zhang (2023) also applied 
control strategies for prosthesis walking on 
stochastically uneven terrain. 

Mechanical actions of muscles are 
necessary for dynamical modeling of the 
human body and can be considered as 
moments applied at the joints. Weiss, Kearney, 
& Hunter (1986) indicated that the moment— 
angle relationships in the hip, knee, and ankle 
joints are similar to that of a non-linear spring. 
Maletsky & Hillberry (2005) designed a 


dynamic knee simulator to reproduce the 
loading and kinematics of the human knee 
during different activities. Lim & Park (2018) 
developed a model for human locomotion with 
a curvy foot connected to a leg by a springy 
segment. Thus, the oscillations of the center of 
mass during walking can be described by the 
mechanics of a simple passive Spring Loaded 
Inverted Pendulum (SLIP). Kim, Lee, & Koo 
(2018) simulated joint reaction forces, active 
moments by muscles, and passive moments by 
connective tissues. They found that, at Chopart 
and Lisfranc joints, passive moments were 
responsible for large portions of the net 
moment. The passive structures and passive 
moments in the midfoot joints provide strength 
and prevent injuries. 

Carrying backpack loads for long distances 
is common in a range of human activities. As a 
result, the influence of backpack carriage on 
physical performance has been investigated to 
establish guidelines for safe load limits. Liu 
(2007) analyzed the effect of backpack load 
position, walking speed, and surface grade on 
the physiological responses of infantry 
soldiers. The results confirmed that positioning 
the backpack mass center as near as possible to 
the body mass center resulted in the lowest 
energy consumption. Alamoudi, Travascio, 
Onar-Thomas, Eltoukhy, & Asfour (2018) 
determined the effect of different carrying 
methods on walking stability using motion 
capture analysis. The results confirmed that 
the lack of stability in the frontal carriage 
forced the body to increase the cadence to 
maintain stability. Additionally, to minimize 
the moment generated by both the upper body 
and the heavy load, participants tended to 
decrease the length of their stride. Walsh, 
Low, & Arkesteijn (2018) studied the effect of 
stable and unstable load carriage on walking 
gait variability, dynamic stability, and muscle 
activity of older adults. The results showed 
that unstable load carriage reduces dynamic 
stability compared to unloaded walking. 

Reviewing the literature reveals that the 
investigation of the joint forces during load 
carrying is important for preventing joint 
injuries. However, to date, no research has 
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been published on the joint forces during the 
working with the motorized backpack 
sprayers. Therefore, the initial objective of the 
present study is to develop a_seven-link 
dynamic model of the operator’s body while 
working with a motorized backpack sprayer. 
The secondary goal of this research is to 
examine the effects of different working 
conditions (engine speed and mass of sprayer) 
on the hip and ankle joint forces. This is a 
continuation of our previous work, where the 
torque—angle relationships in the hip, knee, 
and ankle joints for the backpack sprayer 
operator body were discussed (Karimi 
Avargani, Maleki, Besharati, & Ebrahimi, 
2020). 


Materials and Methods 


In this section, the principle of the Lagrange 
equation for nonholonomic constraints is 
introduced. The total kinetic energy, potential 
energy, generalized forces, and constraints are 
expressed for the proposed model. The non- 
linear equations of motion are formulated 
using Lagrangian equations and solved using 
Maple software (Ver 2015). 


Lagrange Equation for Nonholonomic 
Constraints 


Lagrangian equations have a special place 
in analytical mechanics. They represent 
equations of motion in terms of generalized 
coordinates. A holonomic constraint in the 
system of n generalized coordinates gi can be 
written as Eq. 1. 


& (240-0 Fp rt) =O (1) 
For nonholonomic constraints, Eq. 1 


changes into a differential form in Eq. 2 
(D'Souza & Garg, 1984; Greenwood, 1988): 


> adapt ate TEI, Qos (2) 
v=1 


Where v=/, 2, ..., n and /=/, 2, ..., s are the 


number of coordinates and number of 
constraints, respectively (n>s). Eq. 3 expresses 
Lagrange’s equation for constrained systems. 


d of. OT W<& 
VA, =O). Veh 2 ij (3) 
il 


dt Ga 0g, Og, 


Where 7, V, 4, and Q) are kinetic energy, 
potential energy, Lagrange's coefficients, and 
generalized force related to generalized 
coordinate gi, respectively. These n equations 
have n+s unknowns, namely the n coordinates 
(qv) and the s Lagrange multipliers (4;). The 
additional equations (Eq. 2) are needed for the 
S constraint which is coupled with the q. 
However, as shown in Eq. 4, these are 
considered as differential equations. 


yd +d, =0 l=1, 2, ey S (4) 
v=1 


Dynamic model of operator’s body during 
spraying 


Fig. | shows the proposed dynamical model 
of the operator’s body during working with a 
motorized backpack sprayer. This model 
includes two rigid legs connected to the rigid 
upper body with hinges at the hips. Each leg 
includes the thigh, shank, and foot. The thigh 
and the shank are connected at the knee joint 
and the foot and the shank are connected at the 
ankle joint. A point mass mp at the hip 
represents the pelvis. The mass of upper body 
mp, leg mi, thigh m;, shank ms, and foot mf are 
considered lumped parameters. Longitudinal 
parameters cp, cp ci, Cs, Cr, 1, Ls, and I; are 
distances from the hip joint to the Center of 
Mass (CoM) of upper body, from the ankle 
joint to CoM of foot, from the hip joint to 
CoM of stance leg, from knee joint to CoM of 
shank, from the hip joint to CoM of thigh, 
stance leg length, shank length, and thigh 
length, respectively. The motorized backpack 
sprayer is located at position (Xsp, ysp) from the 
CoM of the upper body. 
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Fig.1. The proposed dynamical model of the operator’s body while working with a motorized 
backpack sprayer 


The total kinetic energy, potential energy, 
generalized forces, and constraints are 
expressed using the following assumptions 
(Wisse, Schwab, & van der Helm, 2004): 

e The joints have no damping or friction. 

e Bones do not suffer from flexible 
deformation. 

e A kinematic coupling has been used in the 
model to keep the body midway between 
the two legs. 

e The shank of the foot stance is always 
locked and the whole leg can be modeled 
as one rigid stick. 

e There is enough friction between the 
walker and the ground. Thus, a flat foot 
does not deform or slip. 

e The muscle moments acting on the hip, 
knee, and ankle joints are modeled as non- 
linear torsional springs. 

e The sprayer engine is an unbalanced 
force. This can be caused by the 
deposition or erosion of the rotational 
parts. 


Kinetic Energy 
As shown in Fig. 1, it is assumed that the x- 
axis is along the ground while the y-axis is 


vertical to the ground pointing upward. The 
dynamic model for the operator’s body while 
working with a motorized backpack sprayer 


can be described by the generalized 
coordinates (Eq. 5). 
gq =x, Yn» Fy, Oy, 5, O55 Oxy J (5) 


Where xn, yn, 07, 02, 3, 02s, and O27 are the 
horizontal coordinate of hip joint, vertical 
coordinate of hip joint, the angle between 
vertical axis and the stance leg, the swing 
angle between vertical axis and the thigh, the 
angle between vertical axis and the upper 
body, the swing angle between vertical axis, 
and the shank and the angle between 
horizontal axis and foot, respectively. The 
positive direction of all the angles is 
counterclockwise. The total kinetic energy of 
the system is defined as the sum of the kinetic 
energy of the point masses, and_ the 
translational and rotational kinetic energy of 
the motorized backpack sprayer. The 
expression for the kinetic energy is obtained in 
Eq. 6. 
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T=T,+T7,+T, +7, +T,, +T, +T, 
1 .2 
= alm +m, +m, +m, +m,, +m, +m, |x, 

1 re 
tlm, +m, +m, +m, +m,, +m, +m, ]¥, 
t m2? +4[ me? +m2 +m,22)6 

peal a) mc, +ML, +Mel, |, 


+ sly +m,c, + m,, iG +¥e y + Xe |e (6) 


1 214 1 ; we ode 
+5[mci + ml? \6;, - 59s +m, (4, cos6, + ¥, sin@) 
—m,c,8; (%, cos 8, — ¥, sin 8, )+,c,8, (%, cos , — J, sin 8, ) 


# m,,; {[(« + yy, a, +X; Joos A, + [yas -(c + sy ) 5, Jsin ax} 


+m, {te,0, cos(0, —,,)—+, (18, cos 8, +¢,6,, cos 6, |+ yy, (16, sind, +¢,6,, sind, } 
+m, {1,1,0,6,, c08( 0, -0,,)+¢;1,6,6,, sin(,, -8,) +¢/1,0,,6,, sin(8,, -8,,) 

-i, (18, cos, +1,0;, cos, +¢,,, sin®, | 

+5,[18,sin8, +1,6,, sind, c,6,, cos, } 


Where J;, is the mass moment of inertia of 
the sprayer about its principal axis. 


Potential Energy 


The total potential energy of the system is 
given in Eq. 7 and is defined as the sum of the 


gravitational energy of the point masses. 
U =U, +U,+U, +U, +U,, +U, +U, 


=-[ml +m, (1 -c,)+m,1+m,1 t+m,,1 +m,1 |g (1-cos@,) 
+[m,c, +m,l,]g (1-cos0,)—| m,c, +m, (c, +V» ) | (1-084) (7) 
+m, gx,, Sin, +m, gc, (1-cos@,, ) 


+m, g (cos, —1, cos, —1, cosO,, —c; sin Ay ) 
Generalized Forces 


The generalized forces acting on the 
operator’s body while working with a 
motorized backpack sprayer include forces 
resulting from gravity and are unbalanced due 
to erosion of the sprayer engine, muscle 
moment acting on the hip joint, muscle 
moment acting on the knee joint, and muscle 
moment acting on the ankle joint. The 
generalized forces Q gi are given by the 
following equations (Eq. 8-14). 


Q,, =mea" cos(at) (8) 
Q,, =mea sin(at) (9) 
C=, (10) 
Q, =T, (11) 
Q, = mea cos at[(c, + y,,)cos@, + x,, sin 8 | 

. +mea sin at[-(c, + y,,)sin @; + x,, cos @ | in) 
Q,, =1, (13) 


Q,,, =0 (14) 

Where mo, e, and @ are the unbalanced 
mass of the sprayer engine, the unbalanced 
mass eccentricity, and the engine speed, 
respectively. 

The muscle moment acting on the hip joint 
T2 (the swing leg, Eq. 15), muscle moment 
acting on the knee joint 72; (the swing leg, Eq. 
16), and muscle moment acting on the ankle 
joint Tif (stance leg, Eq. 17) are non-linear 
functions of the angle between the two 
segments (Karimi et al., 2020). 
T, =0.0007(0, -8,)° -0.0258(0, -0,)° +0.3236(0, —8,) -1.6792 (15) 
T,, =-0.0011(,, -6,)° —0.0029(6,, —0,) +1.2683 (16) 
T,, = 0.00430} + 0.042967 + 0.50528, + 3.1455 (17) 


Constraints 


A constraint can be expressed by a 
relationship between generalized coordinates 
and time. When foot contact with the ground is 
maintained, first and second constraint 
equations can be represented as Eq. 18 and 19. 
y, -Lcos 8, =0 (18) 
x, +1sin 8, — Xin. =O (19) 

Here, the xXankie is the fixed position of the 
ankle of the stance leg. In the passive 
dynamical model of a walking human, the 
upper body can be considered as an inverted 
pendulum jointed at the hip. Therefore, a 
kinematic coupling has been used in the model 
to keep the upper body between the two legs 
and achieve stable walking (Wisse ef al., 
2004). The equation of the kinematic coupling 
constraint is introduced according to Eq. 20. 
20, —(6, +8) =0 (20) 

The time course of hip, knee, and ankle 
joint angles for the proposed model are given 
in Eqs. 21-23 (Karimi et al., 2020). 

0, — 0, =373.64t* — 762.661 +414.86t -38.89 (21) 
0,, —, =-6357.5t? +4911.3t7 —958.48t+1.24 (22) 


6, =451.84t° —513.31t* +218.5t-27.77 (23) 


According to Eq. 4, the constraint equations 
(Eqs. 18-23) can be considered in differential 
form, as defined in Eqs. 24-29: 


j, +1sin 6, =0 (24) 
x, +1cos66, =0 (25) 
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6, +0, -26, =0 (26) 
6, —O, -1120.92t? +1525.32r-414.86=0 (27) 
6, —0,, —19072.5t” +9822.6t -958.48=0 (28) 


6, -1355.52t° +1026.62t —218.5 =0 (29) 

By applying Lagrange’s equation (Eq. 3) 
and introducing the constant coefficients, 
seven coupled non-linear differential equations 


of motion can be obtained (Eqs. 30-36). 
m2, +17, C086, —n; cos 0,0, +[n, cos, +75 sin 8, 16 ~1, 0080.0, 
—1, sin 6, ,,, —n, sin 8,6; +7, sin 0,0; —[7, sin 8, — 7, cos, ]O; 


(30) 


+1, 8in 0,,0;, — 7, cosO, ,O;, +A, = mew cos(@t) 


m5, +7, sin6G, +n, sin 0,6, —[7, sin O, —7, cos, ]6, +7, sin ,,8,, 
—1, c0s6, 6, , +7, cos8,6" +1, cos 0,0; —[7n, cosO, +7, sin 0, |; (3 1) 
+1, c0s0,,6;, +7, sin, 6; + 2, = mea" sin(wt) 

74, +1, (cos 6.x, +sin 4,3}, )—mg sin 8, + A,lsin@, + A,lcos6, + A, + Ag 
=0.00436F + 0.042967 +0.50526, + 3.1455 (32) 


MoO; — 11; C08 A,X, +7; sin A, 5}, +7, c08(O, — A, ),, +1 sin(0,, -0,\6,, 
+m, sin(0, —,,) 0, +1.C08(8,, -0,)0;, +mgsind, +A, +A, +A, (33) 
=0.0007(8, - 0,)’ -0.0258(0, -0,)° +0.3236(0, —O,)— 1.6792 
ns, + (11,c080, +1, sin 9,)x, —(7n,sin8, ~7, c080,)5), — ng sind, 
+7,g CoS 0, — 2A, —A, =mew cosa| (c, +y,,)cosO, +X, sin@, | (34) 
—mew sinat| (c, + y,,)sin 8, —X,, COS 6, | 
Mah, Mg COS, X,, +19, Sin O, %, +1, 008(, — A), )O, + 7,5 sin(0,, -6,, 6, 
—n,, sin(8, -4,,)@; +,:0;, cos(4,, -6,,) +168 sin0,, —A, = (35) 
=~0.0011(6,, —4,)” -0.0029(6,, — 8) + 1.2683 
cy sin, , x, —cos A, , 5, +], sin(0,, -0,)6, +1, sin(0,, -6,,)6,, 
-1, cos( 4, -0,) 0; -l, cos( A, -6,,)0;, —gcos0,, =0 (36) 
The constant coefficients 7; (i=1, 2, ..., 15) 
are defined in the Appendix. Thirteen coupled 
non-linear differential equations (Eqs. 24-36) 
for xn, yn, 91, 82, 63, O25, O2f, 21, 22, 23, 44, As, and 
26 must be solved to analyze the dynamical 
behavior of the operator’s body during 
working with a motorized backpack sprayer. 


Hip and ankle joint forces in the swing 
phase 


Forces in the human joints are important 
factors in the initiation and progression stages 
of joint diseases. Investigation of the joint 
forces has made it possible to prevent or 
minimize discomfort, fatigue, or risk of 
injuries. In this section, the hip and ankle joint 
forces in the model of the operator’s body 
while working with a motorized backpack 
sprayer are calculated. Considering the upper 
body, Newton’s second law leads to Eqs. 37 


and 38. 
DF. =m, x, t,x, +Imye@” cost) (37) 
pee = m, Vi HN,, Vine ts mew sin(@t) (38) 


The components of the hip joint force then 
become Eqs. 39 and 40. 
F. =m,(%, +¢,cos0,0,) +m, (, +X, sin 6,6, + (c, +y,,)cos6,4,) 


X, b Sp ( 3 9 ) 


+me@ cos(at) 
Fy =m,(§, —¢, sin 6,6,) + m,(¥, +X, COS 6,6, -(c,+ y,,)sin 6,6,) (40) 
Finally, the total hip joint force can be 
calculated using Eq. 41. 


i= (EF 41) 


Considering the full body, Newton’s second 
law leads to Eq. 42. 


+myea sin(ot) + (m,, +m,)g 


rs s 2 
in, TF My dy, +X, + MeO cos(at) (42) 


And the x-directional force of the hip joint 
is expressed in Eq. 43. 
F. =m, (x, +¢,cos6,@;) +m, (%, —1, c086,,6,, —1, c08 0,4, —c, sin, 6, -) 

+m,x, + m,(I—c¢,)cos, + m,(%, —c, cos 4,6, —1, cos 0,6) 

+m, (%, +x, Sin 0, + (c, + ¥,)C08,4,) (43) 


+m,(%, —¢, c080,6,) + mew” cos(at) 


Considering the leg stance, the calculation 
of the moment about the hip joint leads to Eq. 
44. 

F, =me(l ~c,)6, + mec, sin8, — F, Leos, /! sin 6, (44) 

Finally, the total ankle joint force is 
calculated using Eq. 45. 


~ 2 2 
P= [F aie (45) 


Results and Discussion 


The values of the physical parameters 
related to this analysis are listed in Table 1. 
Thirteen coupled non-linear differential 
equations (Eqs. 24-36) are solved in Maple 
software. Fig. 2 illustrates the effects of engine 
speed and mass of the sprayer on the hip and 
ankle joint forces during working with a 
motorized backpack sprayer. 

It can be observed that the ankle force 
fluctuates at the beginning of the swing phase 
(t=[0.1-0.2]). It remains relatively constant in 
the middle of the swing phase t=/0.2-0.4] and 
again fluctuates at the end of the swing phase 
(t=[0.4-0.5]). The first fluctuation occurs 
because of the ‘toe-off. In the ‘toe-off 
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instance, the toe loses contact with the ground. 
The second fluctuation was caused by the 
‘heel-strike’. In this instance, there is an 
impact between the leg swing and the ground 
when the heel of the swinging leg comes into 
contact with the ground. The maximum force 
is exerted on the joints when the operator 


begins to swing forward  (‘toe-off’). 
Furthermore, in the instances of ‘toe-off and 
‘heel-strike’, undesirable variations of ankle 
joint force F;, are noticeable. This means that 
ankle injuries due to vibration exposure are 
more likely. 


Table 1- The values of physical parameters (Huang et al., 2012) 


Parameter Value 
Distance from hip joint to CoM of the upper body, cb 0.262 m 
Distance from the ankle joint to CoM of the foot, cf 0.017 m 
Distance from the knee joint to CoM of shank, cs 0.157 m 
Distance from hip joint to CoM of the thigh, ct 0.192 
Unbalance mass eccentricity, e 5mm 
Gravitational acceleration, g 9.81 ms? 
Mass moments of inertia of sprayer about its principal axis, Js, kg.m? 
Stance leg length, / 0.700 m 
Shank length, /; 0.315 m 
Thigh length, /; 0.385 m 
Unbalance mass of sprayer engine, m, 0.01 gr 
Upper body mass, mp 26.62 kg 
Foot mass, my 2.66 kg 
Hip mass, mn 19.97 kg 
Leg mass, my 11.53 kg 
Shank mass, ms 3.45 kg 
Sprayer mass, msp 10.50 kg 
Thigh mass, ™m; 8.07 kg 
Sprayer position from the CoM of upper body, (Xsp, Ysp) (0.175m, 0.038m) 
Engine speed, 3000 rpm 


From a comparison between Fig. 2 (a) and 
(b) it is revealed that the frequency of the hip 
joint force increases with increased engine 
speed (w). However, considerable effects on 
the frequency of the ankle joint force have not 
been observed. Additionally, with increasing 
the engine speed (@), amplitudes of both the 
hip and ankle joint forces increase. 

From a comparison between Fig. 2 (a) and 
(c) it can be concluded that increasing the 
mass of the sprayer (msp) leads to a small 
reduction in the oscillation of the ankle joint 
force. Moreover, as the mass of the sprayer 


(msp) increases, the magnitude of both the hip 
and ankle joint forces increases as well. 

Fig. 3 exhibits the effect of leg length on 
the hip and ankle joint forces while working 
with a motorized backpack sprayer. It can be 
observed that the frequency of the hip joint 
force increases with decreasing the leg length 
l. This is because the leg stance keeps contact 
with the ground while the leg swings and 
pivots about the constrained hip like a 
pendulum. So, shorter operators are more 
vulnerable to injuries due to vibration 
exposure. 
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c) @=3000 rpm, msyp=12.5 kg, 1=0.7 m 
Fig.2. The hip joint force Fy and ankle joint force Fy for 1=0.7 m and (a) msp=10.5 kg, @=3000 rpm, 
(b) Msp=10.5 kg, ©=6000 rpm, and (c) Msp=12.5 kg, @=3000 rpm 
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c) o=3000 rpm, myp=10.5 kg, 1=0.6 m 
Fig.3. The hip joint force Fy and ankle joint force Fr for msp=10.5 kg, o=3000 rpm, and (a) 1=0.7 m, 
(b) 1=0.7 m, and (c) 1=0.7 m 


Conclusion model for the operator’s body while working 
with a motorized backpack sprayer was 


In this study, a novel assistive dynamical : ‘ 
y y presented. In this model, the coordinate of the 
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sprayer relative to the body, rotational inertia 
of the sprayer, muscle moments acting on 
joints, and a kinematic coupling confining the 
upper body between the two legs were 
considered. The dynamics of the sprayer’s 
operator were described using seven 
generalized coordinates. The non-linear 
equations of motion were obtained using the 
Lagrangian equations. The results obtained 
from the numerical analysis indicated that, at 
the beginning and end of the swing phase, 
ankle injuries due to vibration exposure are 
more probable. Moreover, the maximum force 
is exerted on the joints at the beginning of the 
swing phase. Furthermore, the effects of 
engine speed and mass of the sprayer on the 
hip and ankle joint forces were studied. It was 
found that the larger mass of the sprayer (full 
capacity of the tank) can lead to higher levels 
of joint forces and lower oscillations. The 
frequency of hip and ankle joint forces 
increased with the increase of the engine 
speed. The results of this paper can be used for 
an estimated evaluation of a patient’s 
condition and implant design. Investigation of 


Symbols and Abbreviations 


the effects of anthropometric specifications 
and sprayer position (relative to the body) on 
the hip and ankle joint forces while working 
with a motorized backpack sprayer are 
valuable topics for further studies. 


Key Points 


e Development of a seven-link dynamic 
model of the operator’s body while 
working with a motorized backpack 
sprayer 

e The non-linear differential equations of 
motion are formulated using Lagrangian 
equations and solved in Maple software 

e Study of the effects of engine speed and 
mass of sprayer on the hip and ankle joint 
forces 
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Parameter Description Unite 
Cb Distance from hip joint to center of mass (CoM) of upper body m 
Cf Distance from the ankle joint to CoM of the foot m 
Cs Distance from the knee joint to the CoM of the shank m 
Ct Distance from hip joint to CoM of the thigh m 
e Unbalance mass eccentricity mm 
g Gravitational acceleration ms? 
Isp Mass moments of inertia of sprayer about its principal axis kg.m? 
l Stance leg length m 
ls Shank length m 
Lt Thigh length m 
Mo Unbalance mass of sprayer engine gr 
Mb Upper body mass kg 
me Foot mass kg 
mh Hip mass kg 
mi Leg mass kg 
Ms Shank mass kg 
Msp Sprayer mass kg 
m Thigh mass kg 
(Xsp, Ysp) Sprayer position from the CoM of the upper body (m, m) 
o Engine speed rpm 
Xh Horizontal coordinate of hip joint M 
Yh Vertical coordinate of the hip joint M 
a1 The angle between the vertical axis and the stance leg Radian 
&2 The swing angle between the vertical axis and the thigh Radian 
03 The angle between the vertical axis and the upper body Radian 
O25 The swing angle between the vertical axis and the shank Radian 
Oog The angle between the horizontal axis and the foot Radian 
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Appendix 
The constant coefficients for Eqs. 30-36 are as follows: 
7, =m, +m, +m, +m, +m, +m, +m, 1; = mc, 
7), =M,C, + md, + ml, 1, =™M,C, +m,, (c, + Yay) 
Ns; =M,,X,, Ng = M,C, + ml, 
7, = M,C, 7g = mC; 
My = ml ho =mc, +m]; + ml, 
1 = LN hh, = myl.C, 
2 Pa =mc +m,l 
hs = L +m,c, +m,, (e, + Vy) + | ha ss fos 
Ths = M,l,c, 
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